Universal four-body states in heavy-light mixtures with positive scattering length 
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The number of four-body states known to behave universally is small. This work adds a new 
class of four-body states to this relatively short list. We predict the existence of a universal four- 
body bound state for heavy-light mixtures consisting of three identical heavy fermions and a fourth 
distinguishable lighter particle with mass ratio k > 9.5 and short-range interspecies interaction 
characterized by a positive s-wave scattering length. The structural properties of these universal 
states are discussed and finite-range effects are analyzed. The bound states can be experimentally 
realized and probed utilizing ultracold atom mixtures. 
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Universality plays an important role in nearly all ar- 
eas of physics and allows one to connect phenomena 
governed by vastly different energy and length scales. 
A simple class of universal states consists of two-body 
bound states whose size is much larger than any other 
length scale in the problem. Prominent examples in- 
clude diatomic Feshbach molecules [l|, which are nowa- 
days created routinely in cold atom laboratories around 
the world, and di-mesons such as the charmonium res- 
onance near 3870 MeV @- The former have a binding 
energy of order 10~ 10 eV while the latter have a binding 
energy of order 0.5 x 10 6 eV. Yet, once expressed in terms 
of the two-body s-wave scattering length a s , the binding 
energy can be written, to a very good approximation, as 
i5| R rj — fi 2 /(2/ia 2 ) in both cases; here, /x is the reduced 
mass of the constituents (the two atoms and two mesons, 
respectively) . 

Although the concept of universality has been ex- 
tended successfully to three- and higher-body systems 0- 
22j . the list of examples, particularly for few-body sys- 
tems consisting of more than n = 3 constituents, is 
still comparatively small. Most notably, three- and four- 
body physics has been investigated in the context of Efi- 
mov physics. The three-body Efimov effect i.e., the 
existence of infinitely many geometrically spaced three- 
body bound states, can occur when the s-wave scattering 
length is much larger than the range of the two-body po- 
tential. This at first sight purely academic scenario can 
be realized in cold atom experiments by tuning the s- 
wave scattering length in the vicinity of a Fano- Feshbach 
resonance through application of an external magnetic 
or optical field [l|, [l2|]- In the four-body sector, Efimov 
physics can occur via two different routes, as a true four- 
body Efimov effect [22} or as four-body states universally 
tied to three-body Efimov states [fHalsj]- I n either case, 
the description of the Efimov scenario requires two pa- 
rameters, the s-wave scattering length and a three-body 
parameter [j| . The latter may be determined by the char- 
acteristic length of the underlying two-body potential [23| 
(a scenario that has been suggested to apply to cold atom 
systems) or may depend on true three-body physics (a 



scenario more likely to apply to nuclear systems). 

This Letter reports on a new class of universal four- 
body states, predicted to exist — just as Efimov states — 
in three spatial dimensions that are fully determined by 
the two-body s-wave scattering length a s . As such, they 
are fundamentally different from Efimov states, which 
depend on two parameters. Another crucial distinc- 
tion is that the states considered here have continuous 
scale invariance while Efimov states exhibit discrete scale 
invariance. The universal four-body bound states ex- 
ist in heavy-light mixtures that consist of three identi- 
cal heavy fermions and a fourth distinguishable particle, 
which interacts with the heavy particles through a short- 
range two-body potential with positive s-wave scatter- 
ing length a s . We find that the four-body bound states 
exist for mass ratios k larger than k Ci 4 m 9.5. For effec- 
tively two- or one-dimensional confinement, the universal 
tetramers are expected to be more strongly bound than 
in three spatial dimensions. In fact, universal tetramers 
under quasi-two-dimensional confinement have very re- 
cently been predicted to exist for k > 5 [2l|. Just as the 
three-body bound states for positive a s are connected 
to Efimov states (which exist, in the zero-range limit, 
for k > 13.607) [l7|-|S[II H, the universal four-body 
states predicted here are expected to be connected to 
four-body Efimov states, which have been predicted to 
exist for 13.384 < k < 13.607 [Hj]. We analyze the de- 
pendence of the binding energy on the range of the un- 
derlying two-body interaction potential and interpret our 
findings employing hyperspherical coordinates. The uni- 
versal four-body bound states discussed here are not only 
interesting from the few-body point of view but also have 
important implications for the many-body phase diagram 
of heavy-light mixtures |26r l29j that can be realized with 
cold atoms [30M33l | . electrons |34| and quarks. 

Our starting point is the non-relativistic Hamiltonian 
H in free space for n — 1 identical heavy fermions of mass 
M and a single distinguishable light particle of mass m, 
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where 



V tb (r jn ) = -V exp [-r 2 jn /(2r 2 )] 



(2) 



Here, fj denotes the position vector of the jth particle 
and Tjk the interparticle distance, — — Hc|- The 
interaction between the heavy and light particles is de- 
scribed by the Gaussian potential Vtb with depth Vb and 
range r$. In the following, we are interested in the regime 
where the two-body free-space s-wave scattering length 
a s of Vtb is positive and ro -C a s . Throughout, we ex- 
press lengths in units of a s and energies in units of |-Ef R |> 
where i?f R denotes the relative s-wave energy of the two- 
body system with zero-range interactions (realized when 
r -> 0), E% R = -h 2 /(2fia 2 s ) with /i = Mm/(M + m). 
For a given ro/a s , we adjust the depth Vq such that the 
two-body potential supports a single s-wave bound state. 

To determine the eigenstates and eigenenergies of H, 
we separate off the center-of-mass degrees of freedom and 
expand the relative wave function in terms of explicitly 
correlated Gaussians (35|. To construct basis functions 
with good total relative angular momentum L, projection 
quantum number Ml, and parity II, we employ the global 
vector approach 3(| 37} , which has been used extensively 
to describe nuclear systems. Throughout, we limit our- 
selves to states with Ml — 0. The anti-symmetry of 
the basis functions with respect to the exchange of pairs 
of identical fermions is ensured through the application 
of an anti-symmetrizer. The parameters of the explic- 
itly correlated Gaussian basis functions are optimized 
semi-stochastically. According to the generalized Ritz 
variational principle 35[ , the approach yields variational 
upper bounds for the eigenenergies of the ground and ex- 
cited states. Our stochastic variational calculations for 
n = 3 and 4 reported below employ up to 1000 and 3500 
basis functions, respectively. 

We first consider the n = 3 system with L n = 1~ 
symmetry. Employing zero-range s-wave interactions, 
a universal trimer state has been predicted to exist for 
K c,3 ^ 8.173 [2(|. A second universal trimer state has 
been predicted to be supported for k* 3 > 12.917 [20| . 
Symbols in Fig. [TJ show the relative energy £3 of the 
energetically lowest- lying three-body state with 1 _ sym- 
metry, calculated by the stochastic variational approach 
and normalized by |-Ef R |, as a function of ro/a s for var- 
ious mass ratios (k = 8.25 — 10.5, where K = M/m). As 
expected, the trimer energy becomes more negative with 
increasing k for a fixed ro/a s . Moreover, the trimer ener- 
gies approach the zero-range limit from below, with the 
range dependence becoming larger with increasing k. For 
comparison, the dashed line shows the quantity -E2/|-Ef R 1 
as a function of ro/a s ; here, E2 denotes the relative two- 
body energy. Due to the scaling chosen, the dimer energy 
is independent of the mass ratio. The dependence of the 
dimer energy on ro is smaller than that of the trimer 
energy. Dotted lines in Fig. Q] show three-parameter fits 
to the three-body energies with ro/a s < 0.01 38]. The 
symbols in the inset of Fig. [TJ show the extrapolated zero- 




range energies i?f R of the trimer, scaled by |-Ef R |j as a 



FIG. 1: (Color online) Scaled energies as a function of 
ro/a s . The dashed line shows E2/\E 2 R \. The symbols show 
E 3 /\E$ R \ for k = 8.25 toc = 10.5 (top to bottom), in steps of 
0.25; E3 is determined by the stochastic variational approach. 
Dotted lines show three-parameter fits. Inset: Symbols show 
the extrapolated three-body zero-range energy E% R as a func- 
tion of k. The solid line shows a four-parameter fit. 



function of the mass ratio k. The solid line shows a fit 
of the quantity /|-Ef R | to a fourth-order polynomial. 
Our fit predicts that the trimer becomes unbound with 
respect to the dimer for ~ 8.20, which compares fa- 
vorably with the k c 3 yalue of 8.173 determined for zero- 
range interactions 20]. This good agreement demon- 
strates that the stochastic variational approach employed 
in this work is capable of accurately describing universal 
few-body bound states. It should be noted that non- 
universal trimer states can, at least in principle, exist 
for k > 8.6 [l7l l28l. |39(| . Whether or not non-universal 
states exist depends on the details of the underlying two- 
body potential. For the Gaussian model potential consid- 
ered here, it was shown earlier [4Ct |41| that non- universal 
three-body physics comes into play for mass ratios larger 
than those considered here. 

We now discuss the energetics of the four-body sys- 
tem. Circles and triangles in Fig. [5] show the quantity 
E/±I\E!^\ for, respectively, the energetically lowest-lying 
and second lowest-lying states of the four-body system 
with L n = 1 + symmetry for (a) n — 9.5, (b) k — 9.75 
and (c) k — 10 as a function of ro/a s . The four-body 
energies are obtained by the stochastic variational ap- 
proach. Dotted lines are shown as a guide to the eye. 
For comparison, the dashed lines in Fig.[2]show the quan- 
tity i?2/|£ , f R |, and the crosses and solid lines show the 
quantity E^/lE^ R (the symbols show the stochastic vari- 
ational energies and the solid line shows the fit; see also 
Fig. [I). The ground state energy of the four-body system 
lies below that of the three-body system for small ro/a s . 
For k = 9.5, 9.75 and 10, the ground state energy of the 
four-body system "dives down" around r /a s rs 0.015, 
0.011 and 0.008, respectively. In this regime, the four- 
body state acquires non-universal characteristics. For 
slightly larger ro/a s , the energy of the first excited state 
drops below the energy of the trimer and then "traces" 
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FIG. 2: (Color online) Scaled energies as a function of ro/a 3 
for (a) k = 9.5, (b) k = 9.75, and (c) « = 10. The dashed 
lines show i?2/|Sf R | while the crosses (stochastic variational 
energies) and solid lines (fit) show -B3/|-Ef R | (these energies 
are also shown in Fig. [1]). The circles and triangles show 
Ei/\E2 R \ for the energetically lowest-lying and second lowest- 
lying four-body states, respectively. Dotted lines serve as a 
guide to the eye. 



the three-body energy [42| . We refer to the feature where 
the four-body system acquires a new bound state as 
resonance-like feature. The existence and characteris- 
tics of the resonance-like feature depend on the details 
of the two-body interaction model employed. Away from 
the resonance-like feature, the four-body energy shows 
a very similar range dependence as the three-body en- 
ergy, suggesting that the four-body energy is roughly a 
constant multiple of the three-body energy. Moreover, 
it is clear from Fig. [2] that the ratio between the four- 
and three-body energy increases with increasing mass ra- 
tio [43[ . A precise extrapolation of the four-body energies 
to the zero-range limit is challenging for two reasons: (i) 
Numerical issues limit our calculations to ro/a s > 0.004. 
(ii) The existence of the resonance-like feature prevents 
us to perform unambiguous fits. We estimate that the 
four-body system becomes bound around k Cj 4 = 9.5. 

To provide further evidence that the four-body states 
are — away from the resonance-like feature — universal, we 
analyze the hyperradial density P(R). The hyperradius 
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i?cm) 2 + m(f„-i? cm ) 2 , (3) 



provides a measure of the system size [47j . A small hy- 
perradius implies that all n particles are close together 
while a large hyperradius implies that two or more par- 
ticles are far away from each other. In Eq. <j3j) , R cm 
denotes the center-of-mass vector of the n-body sys- 
tem. The hyperradial density P{R), normalized such 
that J °° P(R)dR = 1, indicates the likelihood of find- 
ing the n-particle system with a given R. We calcu- 
late the hyperradial densities as well as other structural 
properties by sampling the n-particle density obtained 
by the stochastic variational approach via a Metropolis 
walk [H[. 

Figures 0(a) and [2Kb) show the hyperradial densities 
P(R) for n = 3 (« = 8.5 and L n = 1~) and n = 4 
(k = 9.75 and L n = 1 + ) for various ro/a s [44|. For these 
mass ratios, the three- and four-body systems support 
very weakly-bound states. To allow for a direct compar- 
ison, dotted and solid lines in the inset of Fig. 03a) show 
the hyperradial densities for n — 3 (n — 8.5) and n = 4 
(k = 9.75) for ro/a s = 0.004. The hyperradial densi- 
ties for n — 3 with k — 8.5 and n = 4 with k — 9.75 
agree qualitatively. They have a small amplitude for 
R/a s <C 1, peak around R/a s = 2 and fall off exponen- 
tially for R^> a s for all ro/a s considered. For fixed k, the 
hyperradial densities move smoothly "outward" with de- 
creasing ro/a s . Importantly, the hyperradial density has 
vanishingly small amplitude not only when R « ro but 
also for notably larger R values (45|. For the three-body 
system, this is consistent with the hyperradial density 
obtained within the zero-range framework (20l . |46| , con- 
firming that the three-body states considered are fully 
universal, i.e., fully determined by a s . The qualitatively 
similar behavior of the n = 3 and 4 hyperradial densi- 
ties for similarly weakly-bound states provides, combined 
with the energetics, strong evidence that the four-body 
states are also universal. 

The dash-dotted line in the inset of Fig.|3Ka) shows the 
hyperradial density of the three-body system for k — 9.75 
and r /a s — 0.004. The three-body system is more 
tightly bound than the four-body system with the same k 
and ro/a s (solid line). In a naive picture, one may imag- 
ine that the four-body system is comprised of a trimer 
with a fourth atom loosely attached to the trimer. Struc- 
tures like this have been observed for the excited tetramer 
state attached to the Efimov trimer state comprised of 
three identical bosons 0, Our analysis of the pair 
distribution functions and radial densities indicates that 
the situation for the tetramers considered here is differ- 
ent. The structural properties of the tetramer of a given 
K loosely resemble those of the trimer with smaller k but 
comparable binding energy. Solid and dotted lines in the 
inset of Fig. EKb) show the radial density n ra d(r J ), nor- 
malized such that 47T J Q n Ta d(rj)rjdrj = 1, for the heavy 
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FIG. 3: (Color online) Structural properties. (a) Solid, 
dashed, dotted and dash-dotted lines show the hyperradial 
density P(R) of the energetically lowest-lying state for n = 3, 
re = 8.5 and r /a s = 0.004, 0.005, 0.006 and 0.015, respec- 
tively. Inset: Dotted, dash-dotted and solid lines show the 
hyperradial density for n = 3 (re = 8.5 and ro/a s = 0.004), 
n = 3 (re = 9.75 and r /a s = 0.004), and n = 4 (re = 9.75 
and ro/a s — 0.004), respectively, (b) Solid, dashed, dotted 
and dash-dotted lines show P(R) for the four-body system 
with re = 9.75 and r /a s = 0.004, 0.005, 0.006 and 0.015, 
respectively. For ro/a a — 0.004,0.005 and 0.006, the energet- 
ically lowest-lying state is considered. For a s /ro = 0.015, the 
energetically second- lowest lying state is considered. Inset: 
Solid and dotted lines show the radial density n ra d(f j) for the 
heavy and light particles, respectively, of the n = 4 system 
with k = 9.75, ro/a s — 0.004 and 1 + symmetry. 

and light particles of the n = 4 system; the position vec- 
tor fj, j = 1, • • • , n, is measured with respect to R cm 
and Tj — \rj\. For large Tj, the radial densities of the 



heavy particles and the light particle are nearly indistin- 
guishable. For small rj, ?i ra d goes to zero for the heavy 
particles but has an appreciable amplitude for the light 
particle, suggesting that the light particle is "shared" 
among the heavy particles. 

In summary, we analyzed heavy-light mixtures in three 
spatial dimensions, where the heavy-light pairs interact 
through short-range potentials with positive s-wave scat- 
tering lengths. Despite the Pauli exclusion principle, 
which acts as an effective repulsion between the identical 
heavy fcrmions, the four-body system consisting of three 
heavy particles and one light particle supports a univer- 
sal bound state if the mass ratio between the heavy and 
light particles is larger than about 9.5. The light par- 
ticle acts as a mediator that "glues" the four-body sys- 
tem together, just as electrons in or H 2 glue together 
the protons by way of the exchange interaction [49[ . Al- 
though the three-body energy shows a fairly strong de- 
pendence on ro, we found that the ratio E4/E3 is, away 
from the resonance-like feature, roughly constant for a 
fairly wide range of ro/a s values, suggesting that the 
universal four-body states can be observed in cold atom 
experiments with current technologies. The existence of 
universal tetramer states opens the possibility to search 
for novel tetramer phases in many-body systems, promis- 
ing a rich phase diagram of heavy-light mixtures on the 
positive scattering length side. In the future, it will be 
interesting to investigate how the universal four-body 
states discussed here are affected by non-universal three- 
and four-body states and how these states are connected 
to Efimov tetramers that have been predicted to exist for 
13.384 < k < 13.607 
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